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We  discuss  a  recently  developed  technique  to  recover  nearly  diffraction-limited  images  of  supergtant  stars  such 
as  o  Ortoms  (Uclctgcusc)  front  speckle  interferometry  data.  This  method  relies  on  the  digital  identification  and 
coaddition  of  the  brightest  individual  speckles  within  a  targe  number  of  shnrt-eaposure  speckle  photographs 
The  resulting  average  speckle  may  be  thought  of  as  the  convolution  of  a  point  source  speckle  profile  with  the 
actual  object  intensity  pattern.  Ily  making  use  of  this  point  we  have  deiived  angular  diameters  and  limb 
darkening  coefficients  in  addition  to  finding  evidence  of  possible  surface  structure  on  the  star.  The  limitations 
of  this  technique  have  been  determined  empirically  and  they  are  discussed  in  this  paper. 
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introduction 

The  application  of  stellar  speckle  interferometry  to 
astronomical  purposes  was  dotte  by  Labeyrte1  and 
Gezari  <•/  <»/ . 1  The  basic  principle  used  In  that  method 
consists  of  measuring  the  mean-squared  modulus  of 
short-exposure  astronomical  photographs.  This  is 
generally  accomplished  by  summing  the  infinity  diffrac¬ 
tion  patterns  of  these  photographs  optically  produced 
using  a  laser.  A  variety  of  theoretical  analyses  of  this 
method1”'  have  shown  that  this  result  contains  recover- 
ible  Information  on  scales  out  to  the  telescope  diffrac¬ 
tion  limit.  However,  since  the  mean-squared  modulus 
is  produced,  phase  information  is  lost  and  a  diffractlon- 
itimted  image  cannot  be  reconstructed. 

The  information  needed  to  completely  define  a  dif¬ 
fraction- limited  image  is  present  within  the  speckle 
data.  Methods  have  been  proposed  and  simulated  to 
extract  this  information. 10  In  such  techniques  the  phase 
information  is  obtained  from  the  statistical  autocorrela¬ 
tion  of  speckle  photographs  In  a  manner  similar  to  that 
used  to  extract  the  modulus  information.  However,  for 
objects  with  sizes  near  the  telescope  resolution  limit 
wo  lave  developed  an  empirical  technique  for  recon¬ 
structing  nearly  diffraction-limited  images."  This 
method  has  been  used  to  study  the  possibility  of  surface 


structure  on  Lite  bright  star  a  Orionis  (Betelgeuse). 
Since  the  detailed  procedures  and  results  are  presented 
In  that  paper  we  will  present  only  a  summary  here  con¬ 
centrating  on  the  limitations  and  subsequent  results. 

EXPERIMENTAL  METHOD 

A  single  speckle  pattern  is  composed  of  individual 
fringes  or  speckles  which  for  a  point  source  object 
have  scales  near  that  of  the  telescope  Airy  disk.5  In 
Fig.  1  speckle  exposures  obtained  on  the  Kilt  Peak 
National  Observatory  4  m  telescope  for  three  different 
types  of  objects  arc  shown.  As  stated  above,  the  In¬ 
dividual  speckles  for  the  point  source  object  y  Orionis 
have  scales  near  the  Airy  disk  size  for  the  4  m  tele¬ 
scope.  Since  a  more  complex  object,  such  as  a  binary 
star  system,  exhibits  a  speckle  pattern  which  is  the 
convolution  of  the  object  intensity  with  a  point  source 
speckle  pattern,  the  speckles  In  the  Image  of  the  binary 
star  a  Auriga  are  each  doubled.  Consequently,  IX  a 
point  source  speckle  pattern  can  be  observed  close 
enough  to  a  program  object  so  that  the  light  from  both 
objects  passes  through  the  same  column  of  turbulent 
atmosphere  (l.c.,  within  the  Isoplanaltc  patch)  then  the 
speckle  pattern  of  the  point  source  object  may  serve  as 
the  Instantaneous  point  spread  function  and  recon- 


KW.  1  •  Speckle  exposuroa  taken  on  tho  Kilt  Peak  National  Observatory  4  m  telescope.  Each  photograph  covers  a  3  arcsoc 
square  area,  (a)  Alpha  Orionis  IBctclgousc),  a  resolved  supergiant  star  with  an  angular  diameter  of  -0.05  arc  sec;  (b)  Gamma 
inionis  tllcllatrix),  a  point  source  star;  (c)  Alpha  Auriga  (Capclla),  a  binary  star  with  an  angular  separation  of  ■-0.0G  arc  sec. 
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strutted  images  may  be  made  using  standard  Fourier 
deconvolution  methods. 11  Dut  there  arc  very  few  suit¬ 
able  point  source  objects  near  enough  Interesting  ob¬ 
jects  (about  2  arc  sec).  ”  Consequently,  the  Lynds 
i-l  nl. 11  technique  relies  on  extracting  point  spread  func¬ 
tion  information  directly  from  the  resolved  object 
speckle  patterns  themselves.  This  is  accomplished  by 
identifying  the  brightest  local  maxima  (i.e.,  the  loca¬ 
tion  of  the  brightest  speckles)  within  the  resolved  ob¬ 
ject  speckle  photographs.  If  this  set  of  positions  Is 
then  assumed  to  represent  the  locations  where  a  point 
source  would  also  exhibit  bright  speckles,  then  it  can 
be  used  as  a  synthetic  "point  spread  function"  to  re¬ 
cover  a  sceing-free  image. 

In  practice  we  use  the  locations  of  bright  speckles  to 
construct  the  profile  of  <m  average  speckle  for  both  a 
point  source  and  resolved  object.  We  then  assume  that 
the  average  object  speckle  o{x,  y )  is  the  convolution  of 
the  average  point  source  speckle  p tv,  y)  with  the  actual 
object  intensity  profile  f(.r,  y),  just  as  the  entire  object 
s|>eckle  pattern  is  the  convolution  of  a  point  source 
S|)ccklc  pattern  with  the  real  object  Intensity  distribu¬ 
tion  i 

o(v,y) -p[x,y) 

Consequently,  by  deriving  o{x,y)  and  pl\',y)  observa- 
uonally,  we  may  recover  l{x,y).  The  key  practical 
limitation  to  this  technique  is  the  ability  to  identify  in¬ 
dividual  speckles.  In  Fig.  1(a),  the  Oetolgeuse  image, 
the  key  limitation  is  evident.  The  fainter  speckles 
overlap  to  such  a  degree  that  It  Is  difficult  to  identify 
individual  ones.  However,  the  brightest  sjiccklcs 
seem  to  be  widely  enough  separated  so  that  they  arc 
individually  identifiable.  For  this  reason  speckle  Idcn- 
tilication  is  restricted  to  the  brightest  speckles  in  each 
photograph. 

A  digital  reduction  scheme  Is  used  to  identify  and  ex¬ 
tract  the  desired  mean  speckle  profiles.  Photographic 
speckle  data  is  obtained  using  the  Kitt  Peak  National 
Observatory  4  m  telescope.11  \yilh  an  Image  scale  of 
2.2  arc  sec/cm,  these  35mm  photographs  have  ex¬ 
posure  times  of  0.01  s.  To  study  the  red  supcrglant 
star  Alpha  Orionis  (Detclgeuse)  40  exposures  of  It  and 
a  comparison  point  source  star,  Gamma  Orionis,  were 
obtained  on  each  of  five  nights  in  March  1974.  Two 
bandpasses  of  100  A  each  were  chosen  for  these  ex^ 
posm  es.  One  at  5180  A  covers  the  TIO  absorption  band 
in  the  spectrum  of  cool  stars  like  a  Orionis.  Light 
from  this  band  represents  a  temperature-sensitive  level 
high  within  the  stellar  atmosphere.  For  comparison 
purposes  a  second  band  at  5100  A  was  used.  The  data 
an  digitized  using  a  PDS  microdcnsltometcr,  with  the 
recorded  photographic  densllics  converted  to  intensities 
using  a  standard  conversion  curve. 

To  identify  the  location  of  the  bright  speckles  each 
speckle  photograph  Is  fit  with  a  two-dlniensional  Gauss¬ 
ian  profile  which  is  then  removed  from  the  data.  This 
serves  lo  delete  a  background  which  consists  of  the  un¬ 
resolved  faint  speckles  and  noise  from  a  variety  of 
sources  such  as  scattered  light  within  (he  system.  To 
unambiguously  Identify  the  center  of  the  bright  speckles, 
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a  Iow-|ias*.digiLil  filter  Is  applied  to  the  data  lo  re¬ 
move  hlgh-frequtfitcy  film  grain  noise.  The  cutoff  of 
this  filler  Is  chosen  to  lie  near  the  telescope  cutoff,  u 
is  then  easy  lo  pick  out  the  bright  speckles  as  the  point, 
of  local  intensity  maxima  brighter  llan  a  prc-clmscn 
threshold.  A  field  of  Impulses,  6(.r,y),  Is  constructed 
which  rep  resents  these  locations. 


fn  order  to  co-add  the  now  located  speckles  in  each 
image  we  consider  each  photograph  to  be  represented 
by  a  simple  model  where  the  image,  0(x,  y),  is  the  sum 
of  discrete  speckles: 

Ok,)-)  =  £/'(*,  y) 

=  £/>(-r,y)»/(x,y);  , 

as  mentioned  previously,  ptv.y)  Is  the  profile  of  a  point 
source  speckle.  This  relationship  can  now  be  re-  "  ' 
written  using  the  impulse  field  6(x,  y)  which  was  con-  ■  ' 
structed  to  represent  the  location  of  the  bright  speckle,. 

0(.v,y)  =  [6(v,y)  •/>(*, y)|  .  /(.v,y).  :  ; 

W c  obtain  our  mean  speckle  profile  by  deconvolving  .  i 
SCv,  y),  which  we  have  derived  from  the  observed  imag,  - 
0(.v,  y).  p(x,  y)  is  obtainable  from  the  point  source  ...  .j 
star  observations.  To  Insure  Dial  systematic  effects  ■;? 
are  minimized  the  point  source  is  chosen  to  He  rela-  1 
tiiely  near  to  the  object  on  the  sky  and  the  exposures  .  5 
are  obtained  tn  time  as  closely  preceding  and  following  >  1 
Ihe  object  exposures  as  possible.  Resulting  mean  ..3 
speckle  closures  averaged  over  all  exposures  obtained! 
during  one  of  the  nights  are  shown  tn  Fig.  2.  *.< 

•I 

Using  these  data  several  types  of  information  may  bei* 
obtained.  Dy  comparing  ihe  images  of  a  Orionis  re-  ...$ 
constructed  in  both  wavelength  regions,  an  indication  of. 4 
surface  structure  Is  present  which  persisted  in  all  flva  f  J 
nights’  observations.  Since  the  results  are  discussed  :ii? 
at  length  by  Lyntls  r(  al. "  they  will  not  be  further  dis-  nj 
cussed  here.  Another  parameter,  useful  for  the  study "tT 
of  conditions  within  stellar  atmospheres,  Is  the  Umb  .  n 
darkening  /(»•),  which  represents  the  change  of  brlght-lm 
ness  on  the  stellar  surface  out  from  the  center  of  the  Jg 
stellar  disk:  It  usually  takes  a  form  ]jj 

f(r)  =  /,[  1  -  v(l  -  cosfl)], 

where  I,  Is  the  Intensity  at  disk  center,  0  Is  the  angle 
between  the  jjoint  at  r  and  disk  center,  and  x  Is  the 
limb  darkening  parameter  taking  values  o*x  « 1.  To 
derive  this  and  the  stellar  radius,  both  assumed  to  be  j 
radially  symmetric  quantities,  the  two-dimensional 
speckle  Images  are  summed  radially  to  produce  the  .ji 
one-dlmenslonal  curves  shown  in  Fig.  2.  It  is  to  be  *1 
noted  that  on  theoretical  grounds,’  although  a  point  f\ 
source  speckle  is  not  expected  to  exactly  resemble  an 
Airy  disk,  the  observed  profile  bears  a  remarkable 
similarity  lo  one,  especially  the  presence  of  secondary  « 
rings.  This  point  source  profile  is  used  to  derive  the  7 
above  parameters  by  convolving  It  with  a  variety  of  |a 
stellar  disk  profiles  varying  the  limb  darkening  and  !■ 
radius  to  provide  the  best  fit  to  the  observed  a  Orlonl, 
profile.  Estimates  of  accuracy  are  acquired  by  de¬ 
riving  these  parameters  using  different  subsets  of  the  , ' 


I244.J 


l.’4J 


Worth'd  cl  al. 


T  Orionis 


t~  .05’  -H 


a  Orionis 
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data.  In  this  manner  very  high  accuracy  (~i  1^)  re¬ 
sults  have  been  obtained  from  a  small  amount  of  data. 

LIMITATIONS 

This  method  is  subject  to  several  sourccsof  systemat¬ 
ic  errors.  The  basic  requirements  of  this  method  may 
be  summarized:  (i)  Discrete  bright  speckles  must  be 
clearly  identifiable  on  individual  exposures,  and  (ti)  a 
comparison  point  source  profile  must  be  obtained  which 
represents  as  closely  as  possible  the  function  which  is 
convolved  with  the  object  to  produce  the  speckles  in  the 
object's  speckle  pattern. 

For  this  technique  to  work  there  must  be  enough 
photons  to  produce  a  speckle  pattern  In  which  the  20  or 
so  brighter  speckles  arc  identifiable.  We  found  that 
speckles  would  not  be  appreciable  smeared  by  the  at¬ 
mospheric  motions  for  exposure  times  as  long  as  0.00 
s.  We  also  found  that  a  spectral  bandpass  of  as  large 
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as  500  A  centered  at  5500  A  could  be  used.  For  larger 
bandpasses  the  higher-order  speckles  near  the  edge 
of  the  pattern  were  badly  elongated,  showing  the  effects 
of  incoherence  over  the  approximately  20  orders  of  in¬ 
terference  present. With  these  limitations  it  was 
found  that  a  speckle  pattern  showing  recognizable 
speckles  could  be  obtained  for  stars  as  faint  as  7th 
magnitude  using  the  Kilt  Peak  photographic  system. 
With  a  more  efficient  TV  recording  system  this  limit 
may  be  extended  to  about  0th  magnitude. 

With  larger  objects  even  the  brightest  speckles  over¬ 
lap  so  that  Individual  speckles  cannot  be  identified. 
Since  there  does  not  exist  a  wide  selection  of  objects 
with  similar  brightness  and  varying  angular  diameter  it 
is  difficult  to  determine  when  this  limitation  becomes 
critical.  However,  for  <r  Orionis,  with  an  angular 
diameter  of  about  twice  the  4  m  Airy  disk  size,  ac¬ 
curate  diameters  which  agreed  well  with  other  deter- 
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TABLE  1.  Derived  angular  diameters  (uniform  disk)  (or  Alpha 
Oi  iuids  obtained  using  differing  roducUon  parameters. _ 


Ddmcicr  (or  a  uniform 
iltslt  (units  of  Q.UUl 
ire  icCumIs) 

Low-pass  filler 
cutoff  frequency 
(uniu  of  airy 
disk  sire) 

Throihold  lor  IdoaUlIrtUo. 
u<  bright  .pee Idea  limit, 
of  Uia  brtght.it  paint  tn 
etch  photograph) 

53  •  4 

Q.G7 

0.65 

52*  3 

1.0 

0.65 

49  .  3 

1.5 

0.65 

53*  4 

1.0 

0.75 

■IS  .  2 

1.0 

0.53 

mi  nations  of  the  diameter  were  obtained.  Recently, 
angular  diameters  and  limb  darkening  coefficients  were 
obtained  for  o  Ceti  (Mira)  which  were  about  three  times 
the  Airy  disk  size.  Clearly  identifiable  speckles  are 
present  for  both  of  these  stars.  No  speckles  are  visible 
in  speckle  photographs  of  the  Jovian  moons  which  have 
angular  diameters  of  about  20  Airy  disks  (-0.6-1 
arc  sec).  However,  for  the  asteroid  Vesta  with  a  diam¬ 
eter  of  about  5-10  Airy  disks,  individual  speckles  may 
be  just  identifiable.  Based  on  these  data  we  conclude 
that  objects  up  to  3-5  Airy  disk  diameters  in  angular 
extent  may  have  images  reconstructed  using  this 
method. 

Implicit  in  our  technique  is  the  assumption  that  each 
speckle  is  an  approximation  to  an  average  speckle. 
However,  by  restricting  the  method  to  bright  speckles 
we  must  consider  the  possibility  that  speckles  sys¬ 
tematically  differ  as  a  function  of  their  brightness.  We 
have  therefore  determined  the  sensitivity  of  our  results 
to  the  method  used  to  locate  bright  speckles  as  well  as 
the  criteria  under  which  bright  speckles  are  defined. 

To  accomplish  this  a  set  of  data  obtained  on  January 
15,  1976  was  reduced  to  determine  angular  diameters 
while  varying  the  parameters  used  to  identify  speckle 
locations.  The  observations  consisted  of  10  speckle 
photographs  each  of  a  Orionis.  These  data  were  re¬ 
duced  as  described  previously. 

We  varied  both  the  cutoff  frequency  for  ihe  low-pass 
/liter  used  to  smooth  that  data  and  the  threshold  in¬ 
tensity  above  which  speckles  were  defined  as  bright. 

The  low-pass  filter  used  In  the  speckle  algorithm  was 
set  at  values  corresponding  to  !,  1,  and  lj  limes  the 
Airy  disk  size  (0.020,  0  030,  and  0.  045  arc  sec,  re¬ 
spectively).  The  threshold  above  which  speckles  were 
called  "bright"  was  altered  between  55%,  65%,  and  75%. 
The  results  of  these  tests  are  listed  In  Table  I. 

It  Is  apparent  that  all  of  the  derived  diameters  agree 
to  within  the  quoted  errors.  However,  the  values  de¬ 
rived  using  a  higher  cutoff  in  the  fitter,  as  well  as  that 
result  which  included  fainter  speckles,  may  be  slightly 
smaller  than  the  others.  This  may  be  attributable  to 
the  effects  of  high-frequency  noise  in  the  photographic 
data.  With  a  filter  that  admits  frequencies  higher  than 
the  telescope  diffraction  cutoff,  substantial  noise  may 
be  passed  causing  some  noise  peaks  to  be  Identified  as 
speckles.  This  would  lower  the  size  of  a  mean  resolved 
object  speckle  proportionately  more  than  the  compari- 
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«u»  point  source  profile  causing  a  decrease  In  the  do 
rived  diameter.  Similarly,  the  lowering  of  the  Inicu^ 
threshold  for  the  Identification  of  bright  speckles  win  j[ 
admit  more  of  the  fainter  speqkles  with  lower  signal.’*, 
to-nolsc.  Consequently,  Including  these  speckles  wogjg 
also  tower  the  derived  diameter  for  the  same  reason.  .1 

It  should  be  noted  that  In  all  cases  where  we  obtain^ 
data  on  a  number  of  different  nights,  the  diameters  do. 
rived  agreed  from  night  to  night  to  wit!. in  the  errors  of 
the  result.  It  Is  concluded  that  our  results  arc  not 
sensitive  to  atmospheric  conditions.  'i 

Wc  have  shown  that  our  method  provides  nearly  dlf.r 
fraction  limited  Images  for  objects  up  to  several  tlme/1 
the  telescope  Airy  disk  size.  Although  there  Is  some  d 
sensitivity  to  high  frequency  noise  in  the  results,  it  jjj 
does  not  appear  to  be  serious.  Current  efforts  are 
underway  to  Improve  the  data  acquisition  system  addlwt 
a  tow-noise  photoelectric  system  which  will  extend  our! 
capabilities  to  fainter  objects.  i’j| 
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